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Abstract. 

Real-time thermal field theory is used to reveal the structure of plasma corrections 
to nuclear reactions. Previous results are recovered in a fashion that clarifies their 
nature, and new extensions are made. Brown and Yaffe have introduced the methods 
of effective quantum field theory into plasma physics They are used here to treat the 
interesting limiting case of dilute but very highly charged particles reacting in a dilute, 
one-component plasma. The highly charged particles are very strongly coupled to this 
background plasma. The effective field theory proves that this mean field solution plus 
the one-loop term dominate; higher loop corrections are negligible even though the 
problem involves strong coupling. Such analytic results for very strong coupling are 
rarely available, and they can serve as benchmarks for testing computer models. 

PACS numbers: 24.10.-i, 52.25.-b 



1. General Formulation 

A nuclear reaction, which we schematically indicate by 1 + 2 — > 3 + 4, takes place over a 
very short distance in comparison with particle separations in a plasma. Hence, it can 
be described by an effective local Hamiltonian density 



The operator K, describes, with interaction strength g, the destruction of the initial 
particles and the creation of the final particles; the operator K) does the reverse. Fermi's 
golden rule presents the rate as 



The angular brackets {•••)& denote the thermal average; Q is the reaction energy release. 

The extension of imaginary time thermodynamic theory to include real time 
behavior was initiated long ago by Schwinger pQ and Keldysh [2j. Using this method as 
a basis, a detailed analysis [3| shows that when the particles entering into the nuclear 
reaction can be treated by Maxwell-Boltzmann statistics, 



W(x, t) = gJC(x, t)+g /C f (x, t) . 
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Figure 1. The C+ portion represents the interactions between the plasma and 
the final reaction particles that appear in (x, x, t\0, 0, Q)Xfif- The C_ part is 
needed for the plasma interactions with the initial reaction particles that enter into 
(0, 0, —i(3Ti\x., x, t)^2 ■ This contour has the real C'_ and purely imaginary C"_ parts. 



with the functional integral definition 

Z c [(f>] = Z- 1 jH [dr b dA] exp J c dsL} 

exp|^dsy(rf 3 y)p(y,s)0(y,s)| . (4) 

All the field variables ijj in the plasma Lagrangian L and plasma charge density p are 
functions of the spatial coordinate y, and the generalized time variable s runs along the 
contour C shown in Fig. 1. The reacting particles have thermal wave lengths A1.2 and, 
with no plasma interactions, they would have number densities nf\. 

The structure of the result (J3J) is easy to understand. The two transformation 
functions (• ■ ■ | ■ • ■) Vcr ^ describe the propagation of the initial and final particles that 
undergo the nuclear reaction. The Vc superscripts indicate that these particles interact 
via their mutual Coulomb forces. The superscripts indicate that these particles also 

produces the 



V 



h 8 



interact with an arbitrary external potential. The operator Zc 
Coulomb interactions between the reacting particles and the background plasma. 

In essentially all cases of interest, one can neglect the real time portions C + and C'_ 
because of the factor exp{iQt/h}: the relevant real time scale is h/Q, a time very much 
shorter than any characteristic plasma time. In many cases of interest, ftr max <C 1, 
where k is the Debye wave number and r max is the turning point radius of the Coulomb 
interaction between the initial particles. Then the rate reduces to j3J EJ El 

<» JV<"» N 1+2 

^=^c lv^lv^]e■ (> 

Here Tc is the nuclear reaction rate for a thermal, Maxwell-Boltzmann distribution of 
the initial (1,2) particles at temperature T but with no plasma background. The rate Tc 
does contain the full effects of the Coulomb forces between the reacting particles. The 
number is the particle number obtained for a free gas grand canonical ensemble 
with chemical potential \x a . The number N a is the particle number of this species a with 
the same chemical potential \i a but now interacting in the plasma. The subscripts 1 + 2 
denote a composite particle of charge [Z\ + Z 2 )e. 
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Figure 2. Ratios of S[i(/ 
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- Zp for the ion sphere model result [short-dashed line] 
and the corrected ion sphere model [long-dashed line] to the exact numerical action. 



2. Method Illustrated By Improving The Ion Sphere Model 



The simplest example has a weakly interacting one-component plasma, g <C 1, where 
g = f3e 2 n/4:7r . The effective field theory of Brown and Yaffe [S] shows that 
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where 



S[x] = / (d 3 r 







(V%(] 



n e 



ie/?x(r) 



1- iefix(r) ~ iZ p e/3 S(r)x(r) 



(6) 
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The normalizing partition function Z is defined by the functional integral whose action 
omits the 8 function term in Eq. (J7J). The tree approximation is given by S[i<j) c \(r)} with 



V 2 0ci(r) 



en 



-/3e<£cl(r) 



Z p eS(r) 



(8) 



This is the familiar Debye-Huckle form, but now placed in a systematic perturbative 
expansion where error can be ascertained. Including the one-loop correction gives 



Det 1 / 2 [-V 2 



K 



N (0) 

p Det 1/2 [- V 2 + k 2 e-^i] 
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We work in the limit where Z p is so large that gZ p S> 1. The point charge Z p e/47rr 
part of 0ci( r ) is large and dominates over a large range. This validates the Salpeter ion 



sphere model which approximates 



1 — exp{— /3e0 c i(r)} ~ 9 (r — r) . The total plasma 



charge in this uniform sphere must cancel the impurity charge, and so Tq = 3gZ p / k 3 . 
The first correction to the leading Salpeter solution can also be computed in analytic 
form except for a numerical integral. Including this correction gives, with C = 0.8498 • • •, 



S[i(f) c 



Z p ^ (3sV 2/3 
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10 ^ lJ ' [- (SgZ p ))j ' 

Brown and Yaffe jH] have shown that the one-loop correction for the background 
plasma with no impurity ions present is given by 

.3 ^ 
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In our limit the term k 2 exp {— /3e<p(r)} in the one-loop determinant can be treated as 
being very slowly varying except when it appears in a final volume integral. Thus, 

Det 1/2 [-V 2 + k 2 ] { k 3 4vr ,1 f 1 



exp -T^T r ° 3 =exp {"3^}' (12) 



Det 1/2 [- V 2 + e-^i] 
This result is physically obvious. The ion of high Z p carves out a hole of radius r in 
the original plasma. The original plasma is unchanged outside this hole. Corrections 
smooth out the sharp boundaries and produce only higher-order terms. The original 
plasma had a vanishing electrostatic potential everywhere, and the potential in the ion 
sphere picture now vanishes outside the sphere of radius r®. Thus the thermodynamic 
potential of the plasma is reduced by the amount that was originally contained within the 
sphere of radius r , and this is exactly what is stated to one-loop order in Eq. (jl2j) . This 
argument carries on to the higher loop terms as well. A term involving n loops carries 
a factor g n . The presence of the impurity modifies this to be Zg n . With g sufficiently 
small, all the higher-order loops make negligible contributions. The corrected impurity 
number N p is hence given by Eq's. (fT2Jl and (jlO}) inserted into Eq. Q. 

The number relation expresses the nuclear rate (j2J) in terms of the tree contribution. 
Including the first correction to the ion sphere result gives 

r = T C exp { 1 (3s) 2 / 3 [(Z, + Z 2 f 3 - ZT - ZT] } 

x exp | V3 C [{Z, + Z 2 f 3 - ZT ~ ZT\ | • (13) 

The first line agrees with the calculation of Salpeter [7j; the second is a new correction. 

The number correction for the number of impurity ions N p placed in the weakly 
coupled background plasma with number N can be used to construct the grand canonical 
partition function Z for the combined system by integrating the generic relation 
N = dlnZ/d/3fi . To simply bring out the main point, we now include only the leading 
terms. Standard thermodynamic relations then lead to the equation of state 

P V= jiV - Z p ^f^- N p \ T . (14) 

Although N p /N may be small, there is a large pressure modification if Z p is large. 
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